Localized state and charge transfer in nitrogen-doped graphene 
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Nitrogen-doped epitaxial graphene grown on SiC(OOOl) was prepared by exposing the surface to 
an atomic nitrogen flux. Using Scanning Tunneling Microscopy (STM) and Spectroscopy (STS), 
supported by Density Functional Theory (DFT) calculations, the simple substitution of carbon 
by nitrogen atoms has been identified as the most common doping configuration. High-resolution 
images reveal a reduction of local charge density on top of the nitrogen atoms, indicating a charge 
transfer to the neighboring carbon atoms. For the first time, local STS spectra clearly evidenced 
the energy levels associated with the chemical doping by nitrogen, localized in the conduction band. 
Various other nitrogen-related defects have been observed. The bias dependence of their topographic 
signatures demonstrates the presence of structural configurations more complex than substitution 
as well as hole- doping. 



Graphene has been proposed as a promising alternative 
to silicon-based electronics for some applications. How- 
ever a reliable control of its electronic properties, for ex- 
ample by chemical doping, is still a challenging taslP^. 
For carbon-based materials, the incorporation of nitro- 
gen in the lattice is a natural choice because of its abil- 
ity to form covalent bonds and to modify the electronic 
properties of sp 2 carbon locally, with minor structural 
perturbations 4 6 . N-doped graphene also offers inter- 
esting prospects for various other applications including 
biosensin^P, field emissiorP, lithium incorporation^^ Q r 
transparent electrodes^. 

The substitution of some carbon atoms by nitrogen 
is expected to give rise to donor states and then to 
n-type doping^^. The synthesis of chemically modi- 
fied graphene has been achieved either by direct growth 
of modified layers 12 15 or by post-growth treatment of 
pristine graphend 7 | 8 | 16 t P^. These studies revealed the 
presence of several atomic configurations for the nitro- 
gen atoms: substitutional ('graphitic'), pyridine-like, or 
pyrrole-like N 8 . 

However, a clear correlation between the synthesis 
methods and the atomic configuration of the chemically 
modified graphene, on the one hand, and between the 
atomic configuration and the electronic properties, on 
the other hand, remains a challenging task. A step in 
this direction has been achieved very recently by Zhao 
et a/P^ through Scanning Tunneling Microscopy (STM) 
and Spectroscopy (STS) investigations of nitrogen-doped 
graphene prepared by chemical vapor deposition on a 
copper substrate with NH3 gas in the feedstock. These 
authors determined the atomic configuration of the nitro- 



gen atoms to be predominantly (90%) a simple substitu- 
tion ("graphitic" nitrogen), with the majority of dopants 
located on the same carbon sublattice of graphene. This 
work left open some key questions relative to the interpre- 
tation of the experimental data. For instance, simulated 
STM images exhibit a depletion above the N atom (see 
also Zheng et a/P) whereas experimental results have no 
such central feature. More importantly, the N-induced 
donor energy level has not been evidenced in their STS 
measurements. 

We present here a STM/STS study of N-doped 
graphene samples on SiC(OOOl) obtained by post- 
synthesis treatment. Our STM images provide clear ev- 
idence of the presence of substitutional nitrogen atoms 
together with more complex structures presenting well- 
defined topographical features. Local spectroscopy, sup- 
ported by simulations based on Density Functional The- 
ory (DFT), reveals for the first time a localized donor 
state in the electronic density of states related to the 
substitutional nitrogen atoms. Based on STM images at 
different biases, we also observed and analysed a change 
in the STM topographic image with respect to the tip- 
graphene distance that reconciles experiment and simu- 
lation. Doping configurations different from substitution 
are also analyzed. 

Our graphene samples were prepared from C- 
terminated n-type 6H-SiC(QQ01 ) w afers, following proce- 
dures available in the literature^^ which lead to mul- 
tilayered graphene with misoriented, decoupled layers^. 
N-doping was achieved by exposure to an atomic nitrogen 
flux produced by a remote radiofrequency plasma source 
fed with N2. The samples were then analyzed using a LT- 
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STM working under UHV conditions (see Supplemental 
Material and reference 11 therein). In our experimental 
set-up, the plasma generator is not in the close vicinity 
of the sample, so that only N radicals (N*) with ther- 
mal energy (and no accelerating voltage) interact with 
graphene. Thus, only the topmost layer is affected by the 
plasma treatment. This post-synthesis doping method 
has the advantage over direct growth of modified layers 
that a well-defined domain could be doped, other part of 
the sample being kept pristine. 




FIG. 1: (Color online) Comparison of STM images (15x15 
nm 2 ) of graphene before (a) and after (b) N* exposure. Inset 
of (a): honeycomb lattice of pristine graphene. Inset of (b): 
2x2 nm 2 image of a nitrogen dopant. Tunneling parameters: 
(a) V s =-0.3 V, 1=15 nA; (b) V s =-0.5 V, 1=500 pA 



Fig. n] presents the graphene surface before (a) and 
after (bj N* exposure. A defect-free honeycomb lattice 
with an interatomic distance of 1.4 A is observed for the 
pristine graphene. After N treatment, many localized 
features appear as bright dots on the images, demon- 
strating the effect of the exposure to nitrogen radicals. 
Moreover, Auger Electron Spectroscopy (AES) confirms 
that the defects are related to nitrogen doping sites (no 
other chemical species (e.g. oxygen) are detected in AES) 
while Raman spectroscopy confirms the presence of de- 
fects sites after surface treatment (see Supplemental Ma- 
terial and reference 10 therein). By visual inspection 
of the STM images, a defect concentration ranging from 
0.5% to 1.1% is deduced, depending on the exposure time 
(0.6% for the figures of the present paper). 

A closer look at Fig. [I] (b) reveals that different dop- 
ing configurations are present. Approximately 75% of 
them display a triangular shape consisting in a bright 
spot (approximately 0.4 nm wide) with a three-fold sym- 
metry (inset of Fig. [I] (b)). These images are very similar 
to the one presented in ref [13| and are assigned to sub- 
stitutional N atoms. No preferential orientation of the 
trigonal pattern is observed and, consequently, both sub- 
lattices of the graphene are affected by the N treatment. 
The graphene honeycomb pattern remains also unaltered 
outside of the vicinity of the defect, demonstrating a very 
local perturbation of the graphene layer. Other typical 
high resolution images of this defect are shown in Fig. [2] 
(a-e) for biases ranging from V§=-0.4 to +0.5 V (see 
Supplemental Material for a more complete overview). 



We also systematically observe that the corrugation at 
the doping site is more pronounced at positive biases, as 
compared to images recorded at negative biases. At this 
stage, the doping of graphene with N atoms can then be 
associated to spatially localized electronic states in the 
conduction band. 
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FIG. 2: (Color online) Topographic images (2.5x2.5 nm 2 ) 
of a substitutional nitrogen atom at Vs=+0.2 V, 1=200 pA 
(a), V s =-0.2 V, 1=100 pA (b), V s =+0.5 V, 1=700 pA (c), 
V s =-0.4 V, 1=100 pA (d), V s =+0.35 V, 1=800 pA (e). In 
the schematic view, the central small dot (blue online) corre- 
sponds to the nitrogen atom, darker and lighter dots (red and 
green online) correspond to the carbon atomic sites around 
which the density of states is higher with decreasing values 
from dark to light (red to green online) as suggested by the 
experimental images, (f)-(g) Simulations for N substitution 
in a 10x10 graphene supercell at Vs=+0.5 V (f) and Vs=-0.5 
V (g) 



Ab initio calculations were performed to obtain the 
local electronic structure and the STM images using a 
Tersoff-Hamann approach. Calculations have been per- 
formed on a (9x9) supercell (0.6% of N atoms) and a 
(10x10) supercell (0.5% of N atoms) with localized basis 
as implemented in the SIESTA package^ (see more de- 
tails in the Supplemental Material and references 12-19 
therein). The simulated patterns (Fig. |2|f-g); see also 
refs [6 13) present a dark (low) spot above the doping 
atom, surrounded by bright dots, corresponding to the 
adjacent C atoms. The central dark spot has been ex- 
plained by charge transfer from the N atom to the neigh- 
boring C atoms which results in a smaller spatial exten- 
sion of the electronic states associated with the N atom 
in the direction perpendicular to the layer compared to 
the one associated with the C atoms forming the C-N 
bonds^. This central 'hole' is not observed in most STM 
images (Fig. [2] and ref 13) but both the weak dependence 
of the pattern shape with the bias voltage and the more 
intense corrugation for positive biases are reproduced by 
simulations. Generally speaking, the limit of our compu- 
tation procedure, besides the intrinsic approximations of 
DFT, are the simplified tip (a 'metallic' V shape atomic 
orbital) and, more important in the present study, the 
small distance between the tip and the atomic layer. For 
numerical reasons, this distance is smaller than 3-4 A in 
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the simulations i.e. close to point contact and underesti- 
mated with respect to the experimental distance. 

Interestingly, the hole at the center of the defect (the 
N atom in chemical substitution) appears in some of the 
experimental images, like Fig. [2[b), i.e. at low negative 
bias (see also Supplemental Material). By varying more 
extensively the experimental conditions, the same pat- 
tern has been observed for Vs=+0.35 V and high current 
(Fig. J2je)). In all these cases (low V s or high current), 
the tip is close to the surface and a better agreement 
with the simulations is expected, and is indeed observed. 
The tip is also closer to the surface for negative bias than 
for positive bias because of the localization of the energy 
state associated with the defect in the conduction band 
(see below). This explains why a plain triangular pattern 
is observed for Vs=+0.5 V whereas a hollow pattern is 
observed at V s =-0.2 V. 
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FIG. 3: (Color online) Comparison of scanning tunneling 
spectra between graphene (black curve) and the simple sub- 
stitution (gray curve, red online), (a): spectra taken with 
the feedback loop active when moving from one spot to an- 
other. Inset: spectra taken with the feedback loop off when 
moving from one spot to another, (b): Topographic image 
of the defect on which the spectra on (a) have been taken. 
Tunneling parameters: Vs=+1.0 V, 1=500 pA. (c) Simulated 
partial-DOS for graphene 9x9 supercell (162 atoms) including 
a single N substitution. The PDOS far away from the nitro- 
gen atom (Cj^k), on the neighboring C atoms (C]_) and on 
the nitrogen atom are represented by the solid, dashed and 
dot-dashed (black, blue and red online) curves, respectively. 



Further insight regarding the local electronic structure 
of the dopant can be gained through STS spectra (Fig. [3] 



(a) with the corresponding image on fig. [3] (b)). The 
spectra display two minima: one at the Fermi energy 
(0 V). a ssociated with phonon-mediated inelastic chan- 
nel d 13 * 25 *, and one corresponding to the Dirac energy at 
-0.5 V. The relationship between the charge carrier con- 
centration n, the Dirac energy Ed and the Fermi Velocity 

(vp) n = 1T (jiv F )'* l ea d s to n = 18. 10 12 electrons per cm 2 
for vp = 10 6 m/s. Assuming 0.6% of nitrogen atoms, 
the charge transfer can be estimated to 0.8 electron per 
dopant atom. From a theoretical viewpoint, our simula- 
tions for a 0.6% doping (Fig.[3](c)) exhibit E D at -0.42 
eV and then a smaller charge transfer (0.55 electron per 
N atoms). These values are larger than those reported 
by Zhao et alf^ and the discrepancy can be explained by 
a number of potential sources of uncertainties: the accu- 
racy in the determination of the nitrogen concentration 
as well as the position of the Dirac point, the presence 
of several types of doping sit es and the uncertainty on 
the value of the Fermi velocit^lIIZl on the experimental 
side or the absence of quasiparticle corrections 33 and 
the regular spacing between the N sites in the supercell 
technique for the DFT calculations. 

The comparison between the dl/dV spectra at the ni- 
trogen doping sites and far from it (Fig. [3] (a)), shows 
that a broad peak centered around +0.5 eV appears in 
the vicinity of the nitrogen atom (dl/dV spectra give 
a quantity proportional to the Local Density of State 
(LDOS)). To the best of our knowledge, this is the first 
experimental determination of the energy level of the lo- 
calized state in N-doped graphene. The simulated Partial 
Density of States (PDOS) located on the N atom, on a C 
atom close to and far from the nitrogen (C\ and C] ou ^ 
atoms respectivily) are displayed on Fig.[3jc). The states 
related to the graphitic nitrogen are clearly obtained in 
the conduction bands but display a double peak structure 
at 0.15 eV and 0.50 eV. The nitrogen-nitrogen interac- 
tions in the periodical structure used in the calculation 
could lead to a splitting of the donor state that would 
not be present for randomly distributed defects 28 . The 
absence of quasiparticle corrections 33 or the reduced 
tunneling current for low bias 2 -^ may also explain why 
the low energy states are not observed. Besides, experi- 
mental STS spectra obtained for a tip above the center 
of the defect probably also probe the PDOS of the 
atoms. Like the N PDOS, the PDOS presents a dou- 
ble peak at 0.15 eV and 0.50 eV but with larger ampli- 
tude of the second one. Moreover a state localized at the 
Ci atoms appears around 1 eV and can be related to a 
second feature in the experimental STS. 

Coming back to the experimental STS spectra, at neg- 
ative bias, the dl/dV signal is found to be lower above 
the nitrogen atom than above the graphene as shown in 
Fig.[3^a). At positive bias both STS spectra have compa- 
rable intensities. To reconcile those measurements with 
the higher corrugation above the dopant site (Figs [l] and 
we have to keep in mind that spectra are measured 
with initial conditions corresponding to the setpoint used 
for the STM image. As a consequence, the tip is located 
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at a larger distance from the atomic plane when a spec- 
trum is measured above the dopant site, compared with 
a STS measurement on the graphene layer and the two 
spectra cannot be quantitatively compared. As it is ex- 
pected that the nitrogen atom lies in the plane of the 
graphene sheet, as demonstrated by DFT simulations 
performed on single^ and bilayer N-doped graphene^, 
the observed corrugation is a purely electronic effect. In 
such cases, the intensity of dl/dV spectra are artificially 
reduced above the N atom due to the higher tip posi- 
tion. In order to overcome this difficulty of interpre- 
tation, we have measured some spectra while scanning 
with the STM feedback loop off, i.e. at a constant tip 
height. In these conditions, STS spectra exhibit a very 
different behavior (inset in Fig. [3] (a)). For negative bias, 
the spectrum recorded above the dopant is slightly larger 
than above the graphene layer, in agreement with what is 
observed on the images. For positive bias, the difference 
is even more spectacular with a LDOS measured at the 
doping site up to five times more intense than the one of 
graphene, strongly confirming that the doping states lie 
in the conduction band. 
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FIG. 4: (Color online) STM images of complex doping sites 
(3x3 nm 2 ). (a), (c) and (e) correspond to three different 
defects measured at a bias voltage of -0.5 V. (b), (d) and (f) 
correspond respectively to the same defects measured at +0.5 
V. 



Besides the substitutional doping, other types of de- 
fects have been frequently observed, examples of which 



are displayed in Fig. [4j Some defects appear higher at 
negative bias (Fig. E] (a) and (c)) than at positive bias 
(Fig. [4] (b) and (d)) contrary to what was observed for 
substitution. Interestingly, a simple substitutional N 
atom is also observed on top of Fig. [4] (a) to (d) close to 
the complex defects. These configurations allow a direct 
comparison of the bias dependence of the simple sub- 
stitutional N (higher at positive bias) and of the other 
defects (higher at negative bias). These new dopant sites 
can then be associated with hole-doping and with a lo- 
calized state lying in the valence band 46 . Furthermore, 
the experimental signature (a triangular symmetry with 
extended oscillations of the wavefunction) is consistent 
with the behavior predicted for a configuration involving 
three N atoms forming a pyridinic configuration^ or with 
a single vacancy^. Other defects, like the one reported 
in Fig. [4] (e) and (f) present a pattern with a rectangu- 
lar symmetry and exhibit a n-type characteristic (higher 
corrugation at positive bias). However, further inves- 
tigations based on STS analysis are necessary to reach 
a definitive conclusion on the atomic configurations of 
these complex doping sites. 

In summary, the N doping of graphene was achieved 
by exposing epitaxial samples to an atomic nitrogen flux. 
The most common (~75%) doping configuration is found 
to be substitutional (graphitic) N atoms. Local STS 
Spectra have provided the first experimental evidence of 
a localized donor state, specific signature of N substi- 
tution. STM topographical images revealed a depletion 
above the N sites characteristic of a charge transfer be- 
tween the N atom and the adjacent neighboring C atoms. 
Finally, other N doping sites with more complex atomic 
structure than a simple substitution were observed, ex- 
hibiting an acceptor-like character. The presented post- 
synthesis treatment opens the door to local tuning of 
electronic properties of graphene which is a prerequisite 
for the development of carbon-based electronics. 
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